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We demonstrate that dynamical probes provide direct means of detecting the topological phase 
transition (TPT) between conventional and topological phases, which would otherwise be difficult to 
access because of loss or heating processes. We propose to avoid such heating by rapidly quenching 
in and out of the short-lived topological phase across the transition that supports gapless excita¬ 
tions. Following the quench, the distribution of excitations in the final conventional phase carries 
signatures of the TPT. We apply this strategy to study the TPT into a Majorana-carrying topolog¬ 
ical phase predicted in one-dimensional spin-orbit-coupled Fermi gases with attractive interactions. 

The resulting spin-resolved momentum distribution, computed by self-consistently solving the time- 
dependent Bogoliubov-de Gennes equations, exhibits Kibble-Zurek scaling and Stiickelberg oscil¬ 
lations characteristic of the TPT. We discuss parameter regimes where the TPT is experimentally 
accessible. 

PACS numbers: 03.75.Ss, 05.30.Rt, 05.30.Fk, 03.65.Vf 


Systems of ultracold atoms provide one of the most ver¬ 
satile platforms for realizing many-body quantum phases 
of matter. In fact, several quantum phases and phase 
transitions such as the superfluid-Mott transition [IHS] 
have been realized in such systems. Yet, many of the 
most interesting phases or phase transitions in such sys¬ 
tems are yet to be observed. One of the most glaring ex¬ 
amples is the elusive antiferromagnetic Neel order [Hi 
in the fermionic Hubbard model, which is believed to be 
a precursor of superconductivity in the model. Another 
example is the recently proposed family of phases based 
on the realization of spin-orbit coupling (SOC) by arti¬ 
ficial gauge fields [ms], which includes topological in¬ 
sulators [TdUIH] . topological superfluids (TSFs) [T7H22], 
and fractional quantum Hall phases [53]. A generic ob¬ 
struction to the observations of many of these phases is 
heating due to spontaneous emission from applied laser 
fields. The heating problem makes it difficult to cool into 
the equilibrium thermal state of many of these topologi¬ 
cal phases. To study these phases, one can also prepare a 
gapped nontopological state and ramp the Hamiltonian 
to drive the system from the nontopological to the topo¬ 
logical state. However, the properties of the short-lived 
topological phase are difficult to probe while it is subject 
to thermal fluctuations. 

In this Letter, we propose a dynamical solution to the 
problem of studying the short-lived topological phase 
by starting the system in its long-lived nontopological 
phase and driving it into the topological phase and back. 
The rapid nature of this process obviates heating; this 
is expected to make our proposal easily implementable 
in experiments. The process involves crossing the quan¬ 
tum phase transition between the phases, which supports 
gapless excitations. Driving through the gapless phase 


transition produces excitations in the gapped phase via 
the Landau-Zener (LZ) transitions [24l [23 with a de¬ 
fect density that demonstrates Kibble-Zurek (KZ) scal¬ 
ing |26H35j . More interestingly, our dip-in-dip-out strat¬ 
egy, where the system is driven through the phase transi¬ 
tion and back, leads to the Stiickelberg interference phe¬ 
nomenon IMIISZ] between the two LZ transitions, which 
in turn results in oscillations of the momentum and en¬ 
ergy distribution of the excitations with the ramp rate. In 
many cases the unique ramp-rate dependence of the ex¬ 
citations’ momentum distributions can be measured via 
standard time-of-flight techniques. This provides an ex¬ 
perimentally viable test for the dynamical fingerprints of 
the topological phase transition (TPT), whose equilib¬ 
rium properties would otherwise be hard to access. 

While this general idea applies to many phase transi¬ 
tions in ultracold bosonic and fermionic systems IMIISHI- 
SOj, we focus on phase transitions whose dynamical prop¬ 
erties are well understood [5011551 HOHiS] . In particu¬ 
lar, we apply this idea to the proposed TSFs [50055] 
in systems of ultracold atoms which host the Majorana 
modes [49053] . Two of the key ingredients [54| for re¬ 
alization of TSFs, namely, controllable Zeeman coupling 
and fermionic Cooper pairing are readily available in cold 
atomic systems. The recent realization of synthetic SOC 
in cold atoms [9HI3] provides the third critical ingredient 
for realizing topological superfluidity thus opening up the 
possibility of observing topological phases in ultracold 
atomic setting. In addition, the challenges of spatial and 
energy-resolved spectroscopy are easily resolved [HESj. 
Despite the advantages of these proposals, the detection 
of TSFs in cold atomic systems is made difficult by the 
low temperature scales involved combined with the heat¬ 
ing associated with SOC. 
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For the one-dimensional (ID) spin-orbit-coupled Fermi 
gases (SOCFGs) studied here, the TFT is accessed by 
raising the Zeeman field past a critical value [TTUI^ [54] . 
Using the self-consistent time-dependent Bogoliubov-de 
Gennes equation (td-BdGE) formalism, we calculate the 
spin-resolved momentum distribution (SRMD) of the 
SOCFGs as it is ramped across the TFT through our 
dip-in-dip-out protocol described earlier. We find that 
the dynamics of the SRMD reflect both Stiickelberg inter¬ 
ference phenomenon and KZ scaling behavior for appro¬ 
priate experimentally accessible ramp rates. We demon¬ 
strate that these oscillations and the scaling behavior per¬ 
sist at finite initial temperature and are robust features 
of the TFT separating the conventional and topological 
phases of the Fermi superfluids (SFs). While a gap clos¬ 
ing is not by itself unique to TSFs, a closing of the gap of 
the nondegenerate Bogoliubov quasiparticles spectrum at 
zero momentum |56j is a yet experimentally unobserved 
smoking-gun signature for a TFT. 

We study ID fermionic atoms with SOC and attractive 
s-wave interactions. The SOG is generated by a pair of 
counterpropagating Raman lasers, with recoil wave vec¬ 
tor kr, energy Er = h^k^/2m, and characteristic time 
scale tr = h/Er^ giving the SOC strength a = Ir^kr/m. 
These lasers couple two hyperfine atomic states repre¬ 
senting the pseudospins a =t, i (for example, |t) = 
1/ = 9/2,mF = -7/2) and |i) = \f = 9/2, = -9/2) 

in atoms [ST]). The transverse Zeeman potential 
strength Dr, set by the Raman coupling strength m , is 
varied in time to drive the TFT. Here we consider vary¬ 
ing Dr linearly from 0 to Dr^ in a time tramp: and back 
in the same time: a piecewise linear ramp protocol of 
duration 2tramp [see blue curve in Fig. [^a)]. Because 
our protocol starts with Raman lasers off (Dr = 0), it 
is straightforward to experimentally realize a long-lived 
conventional SF as the initial state [SS] ; as we will see be¬ 
low, tramp is much less than the system’s lifetime (either 
limited by the spontaneous emission of the Raman lasers 
or inelastic scattering from the Feshbach resonances). 

The system’s Hamiltonian in the Nambu basis 
^k{t) = is H{t) = 

i / (ifc4'fc(t)1''HBdG.fc(t)^'fc(t), where denote the 

annihilation (creation) operators for fermions with mo¬ 
mentum k and spin a. The Bogoliubov-de Gennes (BdG) 
Hamiltonian is miiMI] 

'HBdG,fc(t) = Cfe(t)r^ -I- akr^az + + A(t)Tx, (1) 

where cr and t are vectors of Fauli operators acting 
on spin and particle-hole space, respectively. Here, 
^k{t) = h^k'^/2m — fi{t) combines the kinetic energy and 
the chemical potential /i(t), which is determined self- 
consistently to keep the number of atoms fixed. 

The mean-field pairing potential 

= gir, J (2) 
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FIG. 1 (color online), (a) Time profiles of A(t), and 

Eo{t) for tramp = lOOOtr- The dashed lines denote the times 
whose instantaneous band diagrams are plotted in (b). The 
red dashed lines mark the critical times when TFT happens, 
and the shaded region corresponds to the topological regime. 
Plots are obtained from numer ically solvi ng t he td-BdGE 
[Eq. ( |11[ )] self-consistently [Eqs. ( |12a[ ) and ( |12b| l] with initial 
parameters: Dr(0) = 0, A(0) = 2Er and /r(0) = 0 for SOC 
strength a = 2Er/kr and tramp = lOOOtr-. (b) Quasiparticle 
spectra at different Zeeman potentials Dr. From top to bot¬ 
tom, the energy bands are labeled by E 2 ,k, E^^k, E-i^k, and 
E- 2 ,k- The parameters are as follows: (i) Dr = 0, A = 2Er, 
g = 0, (ii) Dr = 1.56Er, A = 1.93E,., g = —0.02Er, (iii) 
Dr = 2.8Er, A = lAEr, g = —0.14Er, and (iv) Dr = 3.12Er, 
A = 0.91F/r., g = — O.SEr. 


is also self-consistently determined, where (...) denotes 
averaging with respect to the initial thermal distribution. 
The attractive effective ID coupling constant giD < 0 can 
be controlled by Feshbach tuning the three-dimensional 
(3D) scattering length [HMHij . In Eq. Q, we used the 
transformed basis where tpka{t) —f i/kcrit) exp[i^{t)/2], 
giving a real pairing potential: A(t) exp[fz?(t)] —)• A(t). 


The instantaneous quasiparticle excitation spectrum 
of the BdG Hamiltonian [cf. Fig. [^b)] consists of four 
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bands, En,k = sgn(n)e(_i)n where n = ± 1,±2 and 


^±,k 


it) = 




2 1,2 


+A{tr+utr+a^k 


(3) 
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Since HBdG.k respects particle-hole symmetry, the spec¬ 
trum is symmetric around E = 0. As shown in Fig.j^b), 
the instantaneous energy spectrum is gapped for fc 7 ^ 0 ; 
however, for fc = 0 the gap closes when = 

r2R(t)/2 — = 0. Such a gap closing 

without change in the symmetry of the ground state 
(which remains SF for all flu) signifies a TFT [S^IMlfBD] 
between topological [e_,o(i) > 0] and conventional SF 
phases [e_,o(t) < 0]. For flu = 0, the positive and nega¬ 
tive bands are doubly degenerate at fc = 0 ; any nonzero 
riu lifts this degeneracy. 

To study the dynamics around the TPT, we propose 
to prepare conventional SFs [e__o(0 < 0] at nonzero tem¬ 
perature T. We then drive the system through the TPT 
by changing flu according to our ramp protocol with 
riR/ > 2^A^ -I- fij, (where the subscript / denotes the 

quantities at time t = tramp) such that the ramp crosses 
the TPT (cf. Fig.[^. 

We first analytically study the dynamics, considering 
the simple case of slow ramps at T = 0. In this limit, ex¬ 
citations occur near k = 0 and at the transition times t = 
given by the roots of flR(tc) = 2yjAltcY -\- /^(tc)^, 
where the Fermi gas changes from conventional to TSF 
and vice versa. For h^k"^/2m ^ ak, we approximate 


Equation ([^ is a two-parameter driven Hamiltonian |61] 
with instantaneous energy eigenvalues zLEk{t), where 

Ek{t) = \/Eo{t)‘^ -g5(t)2fc2. 

We analyze the dynamics of the TPT using HBdC,k{t), 
where the single-particle state of the system at time t is 
given by 


Mt) = Kit) 


sgn{k)wi^ (t) 




-sgn{k)w^ (t) 

wtit) 


(7) 

with the initial conditions b'^ ( 0 ) = 0 and 6 ^ ( 0 ) = 
1. These two-component vectors are expressed in the 

basis (j)^ with w^{t) = ^^[1 ± E^/t)/Ek(t)]/2. The 

Schrodinger equation for the system then leads to 


ihdtbk(t) = UBdC,k{t)bk{t), 


( 8 ) 


where bk{t) = {b+{t),b^ (t))^. 

We make further analytical progress by ignoring the 
self-consistency condition so that the system can be 
treated as a collection of two-level systems for each 
(/c, —k) pair and use the adiabatic-impulse approximation 
[37l [HlHMj that describes such periodic dynamics accu¬ 
rately for low frequency and/or large amplitude drives. 
Within this approximation, excitations are produced only 
near the critical gap-closing times tc(i, 2 ) when the system 
enters the impulse regime; otherwise, the dynamics occur 
adiabatically in each band and the system accumulates a 
dynamical phase U{tf,ti) = exp[—^ dtEk(t)/h]. In 
the former regime, near the gap-closing times tc(i, 2 )! ex¬ 
citations are produced and the evolution operator is [57] 


E-BdC,k{t) ~ akTzdz 


-f A{t)Tx. (4) 


In this limit, excitations occur only between the Ei^k and 
E-i^k bands [cf. Fig. [^b)]. At fc = 0, the eigenenergies 
are ±Eo(t), where Eq^I) = |-y/A(/pA|r//(t )2 — HR(t)/ 2 | 
with eigenstates 


?o (^) = 






(5a) 

(5b) 


N = \/l-pk[i sin((/?s.fe) - Vz cos((/?s,fe)] - iPy^/Pk, (9) 

where pk = exp (—27r5fc) is the probability of excitation 
formation in each passage through the critical point [241 
[75] with 6k = {ak)'^/{2h\dEQ{t)/dt\t^), and VJs,fc = 7 r/ 4 -|- 
(5fc(ln^fc —l)-|-argF(l—Z(5fc) is the Stokes phase originating 
from the interference of the parts of the system wave 
function in the instantaneous ground and excited states 
at t = tc(i, 2 ) with argr(l —idfc) being the argument of the 
gamma function [TOj. These results give the probability 
of defect formation 

P™ = 4pfe(l-pfe)sin^$st.fe (10) 


where corresponds to positive and negative bands 

[with pseudospin |±) = (|t) ± \\))/V2] and cos9{t) = 
/i(t)/iyA(t)2 -|- /r(t)2. In the subspace of these eigen¬ 
states, the effective low-energy Hamiltonian near k = 0 
is 

'WBdG.fc(i) = a{t)kr]^ + Eo{t)r]^, (6) 

where 5(f) = asin 0 (f), ijx = </o (i)[ 5 o ( 0 ]^ + H.c., 
Vz = ^[/(f)[^[/(f)]1'-^[/(f)[^o (f)]^ and 2r)y = -i[r)z,Vx]- 


at f = 2framp, where $ 31 ,fc = C 2 k + Ps,k is the Stiickelberg 
phase and (/ 2 k = dtEk (f )/h is the dynamical phase 

factor accumulated during passage between the two 
crossings of the gap-closing points |37l HT] |69|. Since 
the excitations occur near /c ~ 0 where the P±i,fc band 
approximately corresponds to pseudospin |±) (along 
the X direction), P™ is directly related to changes 
in the SRMD 6nk± measured along the pseudospin x 
direction. Furthermore, within these approximations, 
|(fPo(i)/rf^k(i, 2 ) = f^R//(2framp), and it can be shown 
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that is a function of ky/t^amp only (see Ref. [51] 
for the derivation). Thus, the integrated change of the 
SRMD 6n± = J dk6nk± displays KZ scaling ~ y/Pamp of 
defect density for a system dynamically evolved through 
the TPT. We now show that these properties persist even 
when the self-consistency conditions for A{t) and /r(t) are 
imposed, as well as at nonzero T (see Fig. [^. 

Change in the SRMD, 


0 0.1 0.2 0.3 0.4 0.5 




FIG. 2 (color online). Change in the SRMD Suk- for spin 
l~) = (It) ~ It)) /%/2 as a function of tramp/tr- and k/kr. For 
large Gamp, the width of the oscillation envelopes scales with 
1 /V^ramp as shown by the red dashed line. Suk- is symmet¬ 
ric with respect to A: = 0; thus, for illustration purposes, we 
only plot Suk- for k > 0. Note that 5nk+ = —Suk-- Inset: 
Integrated change in SRMD 5h_ = J dkSuk- as a function 
of tramp Itr exhibiting oscillations, with the amplitude of the 
oscillations at large tramp scaling like ^tramp, as can be read 
off directly from the y axis. The plots are obtained by numer¬ 
ically solving Eq. (Ill self-consistently [Eqs. (12aI and (12bl] 
with initial conditions fi{0) = 0, A(0) = 2Er, and Dr(0) = 0 
for a temperature AcbT = O.IEf (which is below the critical 
temperature Tc = 0.19Tf |71ll72| i. SOC strength a = 2Er/kr, 
and Dr/ = Z.V 2 Er. 


We solve for the dynamics of the single-particle density 
matrix pf'it) = ('I'j)“(t)'I'^(t)) self-consistently and at 
finite initial temperature, where a, b denote the indices of 
elements in the Nambu basis. The density matrix obeys 
the equation of motion [Eq. 0 ] 

ihdtpk(t) = pk{t)], (II) 

subject to the self-consistency conditions (see Ref. [51] 


for the derivation) 

A{t) = ^ J dkTi{pk{t)Tx), (12a) 

k-it) = J dkTr ipk{t)Ak(t)), (12b) 

where Afc(t) = [‘h?k‘^/2m + akaz) — A{t)Tz. Our sys¬ 
tem begins in the thermal state 

Pk{t) = 'y ' 

n 

E„_k{0)<0 

(13) 

where fn,k = [exp(E„^/j(0)/fcBT)-l-l] ^ is the Fermi func¬ 
tion of the initial Hamiltonian, and Aib is Boltzmann’s 
constant. The wave function Xn,k{i) with its particle-hole 
conjugate Xn.fe(A) = Ty(TyX*_n-k{t) begins as eigenfunc¬ 
tions of the initial Hamiltonian and evolves according to 
i^dtXn,kit) = 'HBdG.fc(A)Xn.fc(A)- Figure a) shows the 
resulting time profiles of the pairing potential obtained 
from solving the td-BdGE (see Ref. [69] for the time de¬ 
pendence of all parameters and remarks on the numerical 
simulation). 

We numerically solved the td-BdGE for the change in 
the SRMD 


d'fT'kdz — Tr I [pfc(2tramp) /^fc(0)] 


^+Tz 


1 ± cr. 


(14) 

Figure [^ shows that Suk- still exhibits Stiickelberg os¬ 
cillations even with inclusion of the self-consistency con¬ 
ditions and at T > 0. Furthermore, for Gamp 
we still see 6nk± ~ fcy'Gamp (see Ref. [55] for an explicit 
demonstration of the scaling), and the integrated change 
in SRMD 6n± = J dkSnk± therefore scales with y/tramp, 
thus, showing the robustness of such interference phe¬ 
nomenon in the present system. We verified that these 
features appear only if Hr/ > 2y^A^ -(- where the 
ramp takes the system through the TPT; thus both the 
KZ scaling and the presence of Stiickelberg oscillations 
mark the TPT. In our calculation, we ignored the effect 
of phase fluctuation as this effect can be suppressed by 
coupling an array of ID SOGFGs [75H75] . 

The parameters used for the plots in Fig. [^ are realistic 
for ID SOCFG experiments. For experiments with 
the Raman laser beams, coupling the |t) = |9/2, —7/2) 
and ID = |9/2,—9/2) states, have laser wavelength 
Xr = 768.86 nm, giving the recoil energy Ep = h x 8.445 
kHz, and time G = H/Er « 20 ps [57]. The single¬ 
body decay time due to photons scattering from the Ra¬ 
man lasers is about 60 ms m, and the lifetime owing 
to three-body recombination is about 200 ms [77]. We 
consider SOGFGs with Fermi energy Ep = Ep. The ID 
Fermi gas criterion is satisfied when E-p < huj^\ for the 
lateral trapping frequency ^>±12^ = 5 x 10"^ Hz, which 
corresponds to characteristic harmonic oscillator length 
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d± = y/hjmxj± « 1345ao, where oq is the Bohr radius; 
the parameters used in the calculation for the plots in 
Fig. i correspond to linear density n « 5 ^ and ID 

interaction strength guy « —O.T^E^Xj. (or 3D scattering 
length 030 ~ —2870ao [SI])- For these values, Fig. 
shows that the Stiickelberg oscillations and KZ scaling 
behavior of the SRMD can be observed within the ex¬ 
perimentally limiting single-body decay time (« 3000tr) 
and thus is feasible experimentally. 

Our dip-in-dip-out protocol is quite general and can 
be gainfully used for observing features related to quan¬ 
tum phase transitions between long-lived and short-lived 
phases of ultracold bosonic and fermionic atoms. In addi¬ 
tion, it provides a route to escaping the heating problem, 
which is one of the major obstacles in measuring proper¬ 
ties of such systems in or near their short-lived phases. 
Moreover, our work also shows that such a protocol ap¬ 
plied to ultracold atom systems, including the one we 
analyzed in detail, may provide us with test beds for 
observation of both KZ scaling P51 - I55] and Stiickelberg 
interference phenomenon |7Mj. 
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Supplemental Material for “Dynamical Detection of Topological Phase Transitions in 

Short-Lived Atomic Systems” 


I. ADIABATIC-IMPULSE APPROXIMATION 


The equation of motion fS9t&fc(t) = 'HBdG,k(,t)bk{t) [Eq. 0], where &fc(t) = (t))^ [Eq. Q] and'HBdG,fc(i) = 

a{t)kr]x Eo{t)r]z [Eq. (|^] with r]^ and rjz being the Pauli matrices acting on the subspace [Eq. ([^], can be 

expressed in form of two-decoupled second-order differential equations as 


J- ih 


Ti9tEo(t) ± EQ{t)dt — 


dta{t) 

a{t) 


ihdt T Eo{t) 




(S-1) 


Assuming no self-consistency, we can use the adiabatic-impulse approximation m to write Eq. ( |S-1| ) as bi,{t) = Vbk{0) 
where the total evolution operator V is decomposed into adiabatic U and impulse N operators. The adiabatic (impulse) 
regime corresponds to the time duration far away from (near) the critical gap-closing time tc(i, 2 )- In matrix form we 








can write down U as 
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Uj = 


3 


0 

o'^Cjk 


J = 1,2,3, 


(S-2) 


where the dynamical phases are given by Cifc = dtEkit)/h, (/ 2 k = dtEkit)/h, and Cafe = dtEkit)/h. 

The impulse operator N can be written as m 


V y/Pk y/l - ) ’ 


(S-3) 


where pk = exp(27r(5fe) is the landau-Zener transition probability |24[ 125) at each critical time, 5k = 
(afc)^/(2/i|di?o(l)/di|t^), ^s,fc = <PS,fe — '^/‘^ and (ps,fc = Tr/d + ^fe(ln(5fe — 1) + argr(l — i5k)- The Stokes phase 
Ps,k increases monotonously from 0 in the adiabatic limit (i5fc —>■ 00 ) to 7r/4 in the diabatic or fast driving limit 
(i5fc —t 0), as seen from the asymptotic argument of the gamma function |70j 


arg r(l - iSk) 


CSk, 5k 1, 

-f-4(ln4-l), 4 » 1, 


(S-4) 


where C « 0.58 is the Euler constant. At the end of the ramp protocol, the total evolution operator becomes 


V = U^NU2NU^=(1^ (S-5) 

\lk Pk ) 

with matrix elements 

Pk = (1 -_pfc)e“*‘^+'= -pfee“*‘^-^ 

7fc = \/(l - + 6“*'^-'=) , (S-6) 

where the phases are given by C+fc = Cifc + C 2 fc + Csfc + 2^s.fc and C-fc = Cifc - C 2 fe + Csfe- The probability of defect 
formation at the end of the ramp protocol (at t = 2tramp) is then given by 

-P™ = \lk\^ = 4pfe(l — Pk) sin^ $st,fc, (S-7) 


where <l>st,fe = C 2 fc + Ps,k is the Stiickelberg phase. Note that in the case of no-self consistency, \dEoit)/dt\ = 
f2Rj,/(2tramp), and consequently 4 is a function of fey^^amp- Since pk and tps.fc are functions of 4, Pk^ is also a 
function of fc-^Aamp- As a result, the defect density displays Kibble-Zurek scaling ^ -^/tramp■ 


II. SELF-CONSISTENCY CONDITION 


The self-consistent chemical potential pit) [Eq. (12b)] is derived from the constraint on the particle density n, i.e., 

/ dkTvfpkit) =n- (S-8) 


Taking the time derivative of Eq. (S-81, i.e., iHdtPkit) = [HBdG,kit), Pkit)], and using the cyclic property of trace, we 
have 


i J dkTrilnBdG,kit), Pkit)]Tz) =0 
i J dfcTr([T^,HBdG,fc(i)]Pfc(i)) = 0 
J dkTriTyPkit)) = 0. 


(S-9) 











Differentiating Eq. (S-9) with respect to time and using the cyclic property of trace, we then obtain 


J dfcTr(Ty['HBdG,fe(f),Pfe(f)]) = 0 

J dkTT{[nBdG,k{t),Ty]pkit)) = 0 
J dkTr [{h^k‘^/2m - p{t) + ak<7z)TxPk{t) - A{t)TzPkit)) = 0. 

Noting that giD / dkTT{pk{t)Tx)/4: = A(t), we then have the self-consistent chemical potential p{t) as 


p{t) = 


giD 

4A(t) 


dfcTr {(t?k^I2m + akaz)TxPk{t) - A{t)T^pk{t)). 


(S-10) 


(S-11) 


III. REMARKS ON THE NUMERICAL SIMULATION 


In the main text, the td-BdGE is given in terms of the single-particle density matrix Pkit). The td-BdGE can also 
be written in terms of the wave function Xn,k{t) = iun,ktit),Un,ki{t),Vn,ki{t),-Vn,ktit)V as 


ikdtXn,k{t) = 'HBdG,/c(f)Xn,/c(f), 


subject to the self-consistency conditions 


m = ^ E / dkT-,k{t\ 


Er,.k{0)<0 


P{t) = 


giD 

4A(t) 


E 






2m 


-E„,fc(0)<0 


where 


Qn,k(t) = {‘^fn,k - 1) '^{\Un,k<T? “ \Vn,ka?), 


(S-12) 

(S-13a) 

(S-13b) 

(S-14a) 

(S-14b) 


with /„ = [exp(£'„ fe(0)/fcBT') + 1]”^ being the Fermi function of the initial Hamiltonian. 

The self-consistent solution of the td-BdGE involves solving a large number of coupled time-dependent differential 
equations (one for each k point). To reduce the number of time-dependent variables, we first calculated the self- 
consistent A(t) and p{t) in the adiabatic regime by solving the time-independent BdG equation. The td-BdGE was 
then solved self-consistently for a small range of states near fc = 0 where excitations occur. Since the ±fc eigenstates 
are related by = a^An^k, we accelerated the computation by focusing on fc > 0. Solving the td-BdGE self- 

consistently with the Zeeman potential varied according the piecewise linear ramp protocol (see blue curve in 


Fig. SI), we obtained a(t), A{t), Eo{t), and p(t) as shown in Fig. 


SI 


IV. SPIN-RESOLVED MOMENTUM DISTRIBUTION 


The change in spin-resolved momentum distribution 6nk- shows Stiickelberg oscillations with the ramp time Gamp 
and for large Gampi duk- scales with ^Gampi as shown in Fig. in the main text. In Fig. S2 we demonstrate the 
scaling more explicitly by plotting Suk- as a function of scaled momentum fc/fc^. i/Gamp/G- 
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- a(t)kr 
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- Eo(t) 

- At(i) 


Supplementary Figure SI. Time profiles of nR(t), a{t), A{t), Eo{t) and for tramp = lOOOtr- The dashed lines denote the 
times whose instantaneous band diagrams are plotted in Fig.[^b) in the main text. The red dashed lines mark the critical times 
when TFT happens, and the shaded region corresponds to the topological regime. Plots are obtained from numerically solving 
the td-BdGE self-consistently with initial parameters: flR(O) = 0, A(0) = 2Er, and /r(0) = 0 for SOC strength a = 2Erlkr 
and tramp — fOOOtr. 



scaled momentum, tramp/tr 



Pi 
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s 

..c 
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Supplementary Figure S2. Change in the SRMD Sn^- for pseudospin |—) = (|t) — |4,)) /\/2 as a function of trampltr and 
klkr^EsaapItr- Note that 5nfc_ is a function of fcyThhnp only for large tramp, as seen from its almost flat nature for small k/kr 
and the width of its oscillation envelopes. The scaling of (5nfc_ can be read off directly from the x axis. Suk- is symmetric with 
respect to fc = 0; thus, for illustration purposes, we only plot Srik- for fc > 0. The plots are obtained by numerically solving 
the td-BdGE self-consistently with initial conditions f2R(0) = 0, A(0) = 2Er and /r(0) = 0 for a temperature /crT = Q.fTlp 
(which is below the critical temperature Tc = O.lQTp | 711 172p. SOC strength a = 2Er/kr, and Hr/ = 3.12i5r. Note that 
5nfe+ = —Srik- due to particle number conservation. 






















